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Abstract - Modularity structures are common in various social and biological networks. However, 
its dynamical origin remains an open question. In this work, we set up a dynamical model 
describing the evolution of a social network. Based on the observations of real social networks, we 
introduced a link-creating/deleting strategy according to the local dynamics in the model. Thus 
the coevolution of dynamics and topology naturally determines the network properties. It is found 
that for a small coupling strength, the networked system cannot reach any synchronization and the 
network topology is homogeneous. Interestingly, when the coupling strength is large enough, the 
networked system spontaneously forms communities with different dynamical states. Meanwhile, 
the network topology becomes heterogeneous with modular structures. It is further shown that 
in a certain parameter regime, both the degree and the community size in the formed network 
follow a power-law distribution, and the networks are found to be assortative. These results are 
consistent with the characteristics of many empirical networks, and are helpful to understand the 
mechanism of formation of modularity in complex networks. 



' Real-world complex networks usually have certain uni- 
versal properties, such as small- world, scale- free, and mod- 
lilarity [IHl] . Scale- free means that the degree of a network 
fiDllows a power-law distribution, and the modularity refers 
• • to the fact that networks typically consist of communities 

. ^ clusters in which the nodes are more highly connected 
l;o each other than that to the rest of the network [T] . In 
;_(the past decade, the investigation on networked systems 

.^iS mainly focused on two aspects. On the one hand, the 
collective behaviors of networked oscillators have been ex- 
tensively studied. Most of these works are carried out on 
static networks, aiming at revealing how network topol- 
ogy affects dynamics [5] . On the other hand, many evolv- 
ing network models have been proposed from the topo- 
logical perspective, aiming at constructing networks with 
power-law distribution of degree [3] and modular structure 
[SIT]- However, from functional perspective, taking biolog- 
ical networks for instance, modularity is generally believed 
to correspond to certain functional groups, where nodes 
within the same group may share similar characteristics 



|5]. This implies that in principle the network topology 
and dynamics are strongly dependent on each other. In 
fact, any formed network structure and dynamical pattern 
are actually the result of coevolution of both network dy- 
namics and topology [S]. For example, in various biological 
and social networked systems, such as the email network 
[10] and the mobile communication network [3], individ- 
uals are more likely to interact with others "similar" to 
themselves [11], which usually leads to networks consist- 
ing of communities driven by shared activities, attributes, 
affiliations, and so on. 

The formation mechanism of modularity and scale-free 
property is crucial to the understanding of structural and 
functional/dynamical properties of complex networks. Re- 
cently increasing attentions have been paid to the adap- 
tive coevolutionary networks [9l fT2l - l24] . For example, Rcf. 
|19| investigated the interaction between link weight and 
dynamical states on networked phase oscillators. How- 
ever, so far, it has not been well understood how modular 
structure and power-law degree distribution concurrently 
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emerge in an evolving network. Motivated by the above 
idea, in the present work, we set up a model of dynam- 
ical network whose nodes are represented by phase oscil- 
lators. Basically, the model is an evolving network which 
grows from a few nodes at the very beginning. The most 
important characteristic of the model is that it describes 
the interplay between the topological structure and the 
dynamics on the network. On the one hand, the node dy- 
namics are coupled to each other according to the network 
topology; on the other hand, the connections among nodes 
can be created or deleted according to the local dynamical 
states. Our particular interest is focused on what types 
of network structures can be formed as a result of the 
coevolution of network dynamics and topology. Mainly, 
our study presented the following results: (i) The com- 
munities, within which the nodes have similar dynamical 
states and the connections arc denser than outside, can 
be naturally formed during the network evolution, (ii) In 
a certain parameter regime, the degree distribution and 
the distribution of community size follow power law. (iii) 
The networks turn out to be assortative, i.e., the nodes 
with high degree tend to connect to other nodes with high 
degree. In addition, the average clustering cocfBcient of 
the network is generally high. These properties exhibited 
by the model concurrently emerge as the network evolves, 
and are consistent with the observations in many realistic 
networks [IHllinilMSHI ■ 

In real networked dynamical system, the individual dy- 
namics is generally complicated and different in principle, 
which is beyond our capability of mathematical treatment 
so far. Therefore, the local dynamics on networks is usu- 
ally simplified as continuous oscillator, or even discrete 
map in theoretical study [Mini HI US] • In our model, 
the node dynamics are represented by phase oscillators, 
which are coupled as in the following dynamical equations: 

N 

dm = t^m +7—2^ 0-mn sin(6'„ - 9„i + 4>mn)- (1) 
'^"^ n=l 

Here, the dynamical state 6 describes the attributes of 
oscillator (node). m,n = 1,2,...,N arc the oscillator 
(node) indices, and 7 is the uniform coupling strength. 
A = {omn} is the adjacency matrix, where = 1 if 
nodes m and n are connected, and amn = otherwise. 
{a;„i} are the intrinsic frequencies of oscillators, and k„i is 
the degree of oscillator m. We noticed that in many so- 
cial networks individuals tend to contact with others with 
similar attributes, i.e., similarity breeds connection pT| . 
Furthermore, if two individuals are in the same environ- 
ment, they arc also more likely to make friends with each 
other. For example, when two individuals study in one 
class, or they work in the same company, they are more 
likely to contact each other [10] . Considering the infiuence 
of environment, in our model we particularly introduced 
an extra phase coupling term (f)mn as 

<Amn = (ipm " ^n) mod (tt). (2) 
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Fig. 1: (Color online) Characterization of the properties of 
networks as a function of uniform coupling strength 7. (a) The 
average cluster size (s) for 4-clique. (b) The relative largest 
community size Rl for 4-cIique. (c) The average clustering 
coefficient (C) . Results are averaged over 100 realizations. The 
parameters are iVo = 5, = 10.0, Ta = 100.0, p = 0.5, and 
the final network size A'^ = 500. 



Here ipm is the average phase of the local order parameter 
which is defined as 

1 ^ 

fm^ '^^ — T ^ ^ ^ran^ ^ • (3) 
n— 1 

In an adaptive networked system, the network topol- 
ogy usually changes according to the dynamical interac- 
tion among nodes during the evolution. Our model is basi- 
cally a growing network which starts from a few seed nodes 
at the very beginning. Particularly, the evolution rule of 
the network incorporates two main manipulations: one is 
the node-adding, and the other is link- adjusting, including 
adding and removing links. It has been shown that in em- 
pirical networks [TU1[H], most new connections are likely 
to build up between one node and its second neighbors, 
i.e., the neighbor's neighbor. Based on this idea, in our 
model we proposed a link- adjusting strategy according to 
the local dynamical states. Depending on the local order 
parameter , an active individual has two options to ad- 
just its link. One option is that if his neighbors do not 
reach a consensus, the individual is free and can link to 
any other individuals. Otherwise, if his neighbors reach 
a consensus, he is "frozen" in his neighborhood and can 
only make new link to his second neighbors. Moreover, ex- 
isting links may be removed from the existing network for 
various reasons [50]. This efi'ect has also been considered 
in our model. 

Specifically, the dynamics and topology of the net- 
worked system coevolve according to the following rules. 

1. At the very beginning, the network consists of Nq iso- 
lated oscillators, whose dynamical states evolve ac- 
cording to Eq. ([T]). The initial states of oscillators 
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are randomly selected from [— 7r,7r], and their intrin- 
sic frequencies are chosen from the distribution of 
g{uj) = 0.75(1 - Lu^), where w G (-1,1). Eq. (P is 
integrated by the fourth order Runge-Kutta method 
with the time step At = 0.01. 

2. At every Ta time (T^ = mAt, m G Af), one new node 
is added into the network by connecting it to an ar- 
bitrarily selected node in the existing network. 

3. At every time (T^ = nAt,n £ Af), one node m, 
randomly selected from the existing network, will be 
activated to adjust his links. Suppose the local or- 
der parameter of node m is r^, which characterizes 



the local coherence around node 



Then with total 



probability rm, the node m first chooses one of its 
nearest neighbors k according to the following prob- 
ability partition function 



select 

n =/K^) 



Lmk{A9,nk) 



(4) 



where dm denotes the set of the nearest neighbors of 
the node to. Lmk{A9mk) = [1 + cos(A0mfc)]/2 is the 
similarity distance [211 [32] • According to the defini- 
tion, the larger the similarity distance is, the closer 
the dynamical states between these two connecting 
oscillators are. Then the node k introduces one of its 
nearest neighbor i to the node to according to prob- 
ability partition function 



add 



n =/(^'*) 



(5) 



In this way, a new link between to and one of its 
second neighbor i is established. Meanwhile, with the 



total probability 1 



the node to connects to an 



existing node in the network according to probability 
partition function f{m,j)j^Q^. 

4. Parallel to the above step, the node to will remove one 
of its existing connections with the total probability p 
according to the following probability partition func- 
tion 

i/Lmk{Ae 

mk J 



n 



(6) 



This implies that the probability removing the link 
between to and k is inversely proportional to their 
similarity distance. 

Following the above network evolution rules, actually 
there are two time scales respectively characterized by 
Ta and Tr in our model. In our simulations, we con- 
sider the situation that adding new node is slower than 
adding/removing links, i.e., Ta is chosen to be larger than 
Tr- In the rare cases where these rules would lead to mul- 
tiple links we do not allow them to be formed. 
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Fig. 2: Typical examples of the final networks: the adjacent 
matrix (top) and the final phase states of oscillators (bottom), 
(a) 7 — 0.2. Approximately, the network is homogeneous and 
the phases of oscillators follow a uniform distribution, (b) 
7 = 6. Communities is formed, and the oscillators within the 
same community have close phase states. The indices of the 
oscillators have been rearranged according to the phase. Other 
parameters are the same as those in Fig. [T] 



According to the above rules, given the parameters, the 
formed network will approximately have a constant aver- 
age degree 

(fc)«2[l + ^(l-p)], (7) 

which is independent on the network size. However, the 
degree distribution can be rather heterogeneous. In or- 
der to investigate the possible communities of the formed 
network, we employed the /c-clique percolation method, 
in which communities can be defined in terms of adja- 
cent cliques [55] . This method has advantage to avoid the 
problem of resolution limit in many other community iden- 
tification methods [11[33]. Specifically, we calculated the 
following quantities to characterize the community struc- 
ture in network: the average cluster size (s) (except the 
largest one), the relative largest cluster size Rl and the 
average clustering coefficient (C) . The average community 
size (s) is defined as 



is) 



(8) 



where s is the size of cluster and Hs is the number of s- 
size cluster. The sums run over all possible values of s but 
the largest cluster. The relative largest cluster size Rl is 
defined as 

= ^, (9) 



where sl is the size of the largest cluster, and N is the 
size of system. In our model, the active node can con- 
nect to one of its second neighbors with a certain prob- 
ability, which generates at least one triangle. Therefore, 
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Fig. 3: (Color online) The topological properties of final networks oi N = fOOOO. (a) The distribution of degree. Note that 
the inset plot is on log-normal coordinate, (b) The distribution of community size for 4-clique, where 7 = 6. (c) The average 
nearest neighbor degree fc„„(fc) as a function of degree fc, where 7 = 6. (d) The average clustering coefficient C{k) as a function 
of degree k, where 7 = 6. The straight line is a guide to the eye. Other parameters are the same as those in Fig. [1] 



we focused on 4-clique and 5-clique in our computation. 
Actually, our numerical experiments show that they give 
qualitatively consistent results. 

Now we report the main results of our numerical simu- 
lations. We first study how the network evolves with the 
increase of the coupling strength. Figure [1] shows vari- 
ous properties of the final networks with respect to the 
coupling strength 7. It is found that when 7 — 0, both 
(s) and Rl are very small. This means that the network 
structure is basically homogeneous and no obvious clusters 
are formed at this stage. However, the network structure 
significantly changes when 7 is large enough. For instance, 
when 7 > 2 , both (s) and Rl are significantly large in- 
dicating that the network structure is characterized by 
limited number of communities. This can be further ver- 
ified by the average clustering coefficient (C) as shown in 
Fig. [IJc). A careful examination of the dynamical states 
on the network shows that when the coupling strength 7 
is small, the system can reach neither synchronization nor 
clustering. In this case, the links are almost randomly 
generated, so the formed network is a random one. There 
are no obvious communities as shown in Fig. [DJa). On 
the contrary, when 7 is sufficiently large, clustering occurs 
in the networked system with the growth of the network 
size. Physically, there are two factors affecting the for- 
mation of dynamical groups. On the one hand, according 
to the network-growing rule, a node is more likely to con- 
nect to its second neighbors with similar states. This will 



gradually generate a core of oscillators which are partially 
synchronized. On the other hand, if two oscillators are in 
two different dynamical groups, their "environment" are 
different, i.e., 0„iti between them is large j34j . Since large 
(pmn is not in favor of synchronization |35j , these two nodes 
have less chance to build a connection between them. In 
addition, once the phase difference between two oscillators 
is large, it is likely that the existing link between them 
will be disconnected, while new link is hardly generated 
between them. Therefore, once a dynamical group, i.e., a 
core of oscillators with similar dynamical states and sim- 
ilar frequencies, is generated, the network-growing mech- 
anism of the present model will enhance the formation of 
dynamical groups, as well as the formation of dense con- 
nections inside the group. Numerically, it is observed that 
after a long time evolution, the oscillators self-organize 
into many communities, both dynamically and topologi- 
cally as shown in Fig. [Hb). In our model, the coupling 
strength 7 stands for the magnitude of strength of interac- 
tion among different nodes in the network. Therefore, the 
above results emphasize one important fact that during 
the evolution of networked system, the network structure 
can be significantly affected by the interaction strength 
among nodes. Conversely, if we want to correctly ana- 
lyze the formation mechanism of modularity in a network, 
we have to pay more attention on the specific dynamical 
processes on it. 

In Fig. 121 we particularly provided two examples to 
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Fig. 4: (Color online) Characterization of tlie influence of parameters on the properties. Average clustering coefficient (C), 
relative largest community size Rl and average community size (s) as a function of Ta/Tr (a), as a function of p (b) and as a 
function of A'^. Results are averaged over 100 realizations. Other parameters are the same as those of figure [T] 



demonstrate how the network structure, such as the degree 
distribution and the community size distribution, varies 
with the uniform coupling strength 7. For small coupling 
strength, e.g., 7 ~ 0.2. it is found that basically there are 
no communities as shown in Fig. [2] (a). This is because 
of the fact that when the coupling strength is small, the 
interactions among nodes are weak. Thus it is difficult to 
form dynamically coherent state on the network. In this 
case, most new links are generated randomly. As a result, 
the network is of random nature and the degree approxi- 
mately follows an exponential distribution as shown in Fig. 
[3] (a). Nevertheless, when the coupling strength becomes 
large enough, e.g., 7 > 2, the formation of communities 
is greatly enhanced by the dynamical interactions among 
nodes on the network. In this case, we observe approx- 
imate power-law distribution of degree as shown in Fig. 
[3] (a) in a certain parameter regime. In Fig. [3] (a), it is 
also seen that the power-law distribution has the "droop 
head" shape, which has been observed in many empirical 
networks [3]. In our study, it is found that the distri- 
bution of community sizes also approximately follows a 
power law as shown in Fig. [3] (b), a phenomenon which 
has been observed in many empirical networks |25[|28j . In 
particular, we computed the average nearest neighbor de- 
gree knn,i = '^ij^j characterize degree-degree 
correlations in the network [26| . Averaging this quantity 
over the whole network gives the averaged nearest neigh- 
bor degree fc„„(fc). As shown in Figl3](c), its distribution 
turns out to be power-law with positive exponent. This 
obviously demonstrates that the formed network is assor- 
tative. Moreover, we also computed the average clustering 
coefficient C[k). It decays with the degree k as shown in 
Fig. [3] (d). This property has been commonly found in 
many empirical networks |26ll27j . 



In this work, we have further investigated how the 
formed network structures depend on the main control pa- 
rameters, i.e., Ta/Tj., p, and the final network size N. The 
results are presented in Fig. HI As shown in Fig. IDJa), the 
parameter Ta/Tr has remarkable effect on the properties of 
the formed networks. The main network quantities, such 
as (s), Rl, and C, all increase with the increase of Ta/Tr- 
In our model, the ratio Ta/Tr represents the interaction 
between two opposite regulations. Ta is the time period 
to randomly build a link, while Tr is the time period to 
"purposely" build a link. The former manipulation is in 
favor of increasing randomness in the network, while the 
latter manipulation helps generate triadic closure inside 
the network. Larger Ta/Tr means that the second factor 
is dominant, and this explains why network quantities (s), 
Rl, and C increase when Ta/Tr becomes larger. In our 
model, the parameter p is the probability deleting links 
with small similarity distance. When p is relatively small, 
this manipulation helps form dynamical groups, where os- 
cillators within the same group have similar dynamical 
states, and this process in turn enhances the formation of 
topological communities in the network. However, with 
further increase of p, according to Eq. ([7]), the network 
gradually becomes too sparse to form any distinct commu- 
nities. Therefore, there exists an optimal parameter range 
0.2 < p < 0.4, where communities with various sizes can 
be formed in the network. In our model, the final network 
size N is proportional to the total integration time. As 
shown in Fig. [3)Jc), all three network quantities (s), Rl, 
and C become approximately stationary when the network 
evolves for a long time. Particularly, it is found that the 
clustering coefficient (C) in the formed networks is gen- 
erally large (near 0.3). This property is consistent with 
empirical observations [^[?7]. 
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To summarize, we have studied a network-growing 
model of phase oscillators, in which the dynamics and 
the network topology interact with each other and con- 
currently evolve. Following simple node-adding and link- 
adding/rcmoving rules, the model exhibits several inter- 
esting behaviors. Within a certain parameter range, the 
dynamical communities and the topological modules can 
spontaneously emerge in the network. It is found that 
in the formed networks, the degree and the community 
size approximately satisfy power-law distributions. Fur- 
thermore, the formed networks show assortative connec- 
tion patterns and exhibit high clustering coefficients. Our 
study also reveals that the interaction strength among 
nodes on network can essentially determine the formation 
of network structures, both dynamically and topologically. 
This is an important point which has been ignored or 
over simplified by many previous network-growing models. 
The findings in this work capture the typical properties of 
many realistic networks. Thus it is helpful for us to fur- 
ther understand the complicated interactions between the 
network topology and the dynamics. 
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